Fermi excitations in a trapped atomic Fermi gas with a molecular Bose condensate 
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We discuss the eflFect of a molecular Bose condensate on the energy of Fermi excitations in a 
trapped two-component atomic Fermi gas. The single-particle Green's functions can be approxi- 
mated by the well-known BCS form, in both the BCS (Cooper pairs) and BEC (Feshbach resonance 
molecules) domains. The composite Bose order parameter A describing bound states of two atoms 
and the Fermi chemical potential ^ are calculated self-consistently. In the BEC regime character- 
ized by < 0, the Fermi quasiparticle energy gap is given by y/i^ + A^, instead of |A| in the BCS 
region, where ^ > 0. This shows up in the characteristic energy of atoms from dissociated molecules. 
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There has been increasing interest in the BCS-BEC 
crossover in a two-component trapped atomic Fermi 
gas0, 0, This recent theoretical work has extended 
the classic work in superconductors00,0,|3| to include 
a Feshbach resonance]^ . This allows one to both in- 
crease the attractive interaction|^ |^ by working above 
the resonance as well as to include the effect of long- 
lived dimer molecules just below the resonance. Several 
groups'!^ 0, ^2 have presented evidence for a molec- 
ular Bose condensate in a two-component Fermi gas in 
the region af^ > 0, where af^ is the two-body s-wave 
scattering length. Recent theories fll 0, 1^ clearly show 
that a unified description of both the BEC and BCS lim- 
its can be given, in terms of a composite order parameter 
which involves a Cooper pair condensate and a molecular 
Bose condensate associated with the Feshbach resonance. 
They are two limits of a superfluid Fermi gas, both aris- 
ing from a Bose condensate of bound states. The only 
difference lies in the origin of these bound states. 

This similarity has been noted in the study of the 
collective modes of this condensate in the BCS and 
BEC limits [T^ IT^ as well as in the crossover region^. 
In the present letter, we emphasize this similarity by 
considering the single-particle Fermi quasiparticle spec- 
trum in the BCS-BEC crossover region with a Fesh- 
bach resonance. We point out that the well-known BCS- 
Bogoliubov excitation spectrum is a good approximation 
in both the BCS and BEC limitsi2j at all temperatures 
in the superfluid phase. In particular, the Fermi excita- 
tions exhibit an energy gap due to their coupling to a 
molecular Bose condensate, the analogue of what hap- 
pens in the classic BCS limit. Our work emphasizes the 
importance of measuring|l5j the energy and momentum 
of single-particle excitations in trapped superfluid Fermi 
gases as a way of studying the effect of the molecular 
Bose condensate. We do not address possible corrections 
very close to the Feshbach resonance in the unitarity limit 
(fcpaf ^ ±oo). 

We consider a gas of Fermi atoms composed of two 
atomic hyperfine states (labeled by a =t,i), coupled 



to a molecular two-particle bound state. The coupled 
fermion-boson model Hamiltonian is given bv[lLl3.l9lllq| 



H 



pa q 

^ X! 4+q/2T'^-p+q/2i'^-P'+q/2iV+q/2T 

p,p',q 

9r^[blc^P+q/2iCp+q/2] +I1.C.]. 
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For simplicity, we discuss the case of a uniform two- 
component Fermi gas, although results in the figures 
take into account an isotropic harmonic trap poten- 
tial. In Eq. Q, a Fermi atom and a Boson associ- 
ated with the Feshbach resonance are, respectively, de- 
scribed by the destruction operators Cpo- and 5q. The 
kinetic energy of a bare Fermi atom is = /2m and 
^Bq + 21/ = /2M + 2v is the excitation spectrum of 
the bare 6-molecular bosons. The lowest energy of the 
6-bosons (2v) is the threshold energy of the Feshbach 
resonance. The last term in Eq. is the Feshbach 
resonance with a coupling constant gr, which describes 
how a 6-molecule can dissociate into two Fermi atoms 
and how two Fermi atoms can form one 6-boson. Eq. 
also includes the usual attractive interaction — [/ (< 0) 
arising from nonresonant processes. 

Since a 6-Bose molecule consists of a bound state of 
two Fermi atoms, the boson mass is M = 2m and the 
conservation of the total number of particles TV imposes 
the relation 



TV : 



E(4 



po-Cpo-/ 
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We incorporate this crucial constraint into the model 
Hamiltonian in Eq. JQ) using a single chemical poten- 
tial and work with 1-1 = H ~ fjiN . This leads to a shift in 



the energies, £p ^p 
esq + 21^ - 2/i. 



/i and esq + 2v ^ 
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The mean-field approximation (MFA) for the coupled 
Fermi-Bosc Hamiltonian reduces tof 



/[FA)fo: 



per p 

+ ^^Bcibl^bq + g,^[(l),nCpiC-p^ +h.c.],{3) 

q#0 p 

where the Cooper pair order parameter is A = 
[/ ^p(c_p|Cp|) and the molecular condensate is de- 
scribed by (pm = {bq=o). Thus the MFA broken sym- 
metry state is described by 

per q^O 

- Jlt^'^Pi'^-pT+M, (4) 
p 

where the composite order parameter is given by A = 
A 

~ 9r4'm- The molecular Bose condensate (f)m and the 
Cooper pair order parameter A are strongly hybridized. 
One can also show thatd H H H 



which leads to a renormalized pairing interaction UcS 

^2 



(5) 



piCpt), C/cff = [/ + 



5r 



2v ~ 2fi 



(6) 



In recent discussions on the BCS-BEC crossover (for 
example, Refs. [11 [13), the Feshbach resonance is of- 
ten only included to the extent that the order parameter 
A in Eq. ^ involves the correct two-body s-wave scat- 
tering length given by —Airaf^h'^/'m. = U + gf./2v. The 
two-body Feshbach resonance at 2;^ = occurs at the 
resonance value Bq of the magnetic field. Such calcula- 
tions are identical to the original theories of the BCS- 
BEC crossover in superconductors |3, [E [3- explicit 
account is taken of the formation of a ^-molecule con- 
densate described by as done in the above MFA cal- 
culation, which leads to the effective pairing interaction 
in Eq. ®Q,I1I3- 

Using the simple pairing MFA summarized by Eqs. 
Q-® gives the usual BCS-Gor'kov expressions for the 
diagonal and off-diagonal single-particle Green's func- 
tions 



Gii(p,w) 



^ + ^P 



G'i2(p,w) 



(7) 



with the BCS-Bogoliubov excitation energy given by 



Ep = V(£p-/^)^ 



lAP. 



(8) 



Using (712 (p, to calculate A, one finds that A satisfies 
the "gap equation" with the pairing interaction J/cff, 



(9) 



As usual, a cutoff is introduced in the momentum 
summation |3J. At this MFA level, the total number 
of atoms at T = is 



(10) 



where the number of Fermi atoms is given [using 
Gii(p,iu;)] by the well-known expression 



1 



•per/ 



p,cr 



^[l-^tanh-/3i?p. 



(11) 



The number of Bose-condensed 6-molecules is — 
|0mP, where 4>rn is determined by A and ^. 

We briefly sketch how one proceeds in the case of 
an isotropic harmonic trap(2Cl|. We expand the fermion 
quantum field operator ^'cr(r) and boson operator $(r) 
in terms of the single-particle eigenfunctions /„;„i(r) = 
Uni{f')Yi„i{0 , (j)) of the harmonic potential, with energy 
= {2n + l + 3/2) awo- The resulting MFA Hamilto- 
nian, which corresponds to Eq. Q for a uniform gas, is 
given by [13 

nlm,(7 nlm 

- E PL'[cii„.A'i-mi + h.c.]. (12) 
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Here c„;„j {bi„in) is the creation operator 
of an atom (molecule) in the state fnim{Y) 
Unimi'^))- ^nn' ^rises from the pair potential, 
Kn' = r'^drA{r)u„i{r)u„'i{r). Eq. ^ can be 
diagonalized by the usual BCS-Bogoliubov transforma- 
tion. One can calculate A(r) ~ L/(\E'|(r)^'f (r)) as well 
as A^F in Eq. by expanding the fermion operators 

in terms of the BCS-Bogoliubov excitation operators. 
Using the generahzation of Eq. Q, 



-A(r) 



\-%^ + \Muy+2v^2^ly,r.{r)^0, (13) 



we can determine 0m (r) from A(r), and hence find the 
composite order parameter A(r) in as well as N^. 

The values of /x and A are determined by the self- 
consistent solutions of the MFA gap equation © and 
the number equation given by Eqs. (|10|l and 111|) . This 
simple "pairing approximation" for the single-particle ex- 
citations is expected to give a quantitative description 
at T = (where tanhi^Sp 1) since fluctuations 
are small and all the 6-molecules are Bose-condensed. 
This T = limit was first studied by Leggett|3| in the 
absence of a Feshbach resonance. For T approaching 
Tc, however, the fiuctuations associated with exciting 
^-molecules out of the condensate and coupling to the 
particle-particle (Cooper-pair) channel become increas- 
ingly important 1, 3]. These fiuctuations (rather than the 
breaking up of two-particle bound states) determine 
and the region above, as first discussed by Nozieres and 
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FIG. 1: The number of particles as a function of the threshold 
energy 2v in a trapped superfluid Fermi gas at T = 0. The 
Fermi energy ef = 31.5?i(jJo is for a free gas with = 10912 
atoms. Figs. 1-3 are for U — 0.52eF and gr = 0.2£f. The 
inset shows UeS in a uniform Fermi gas, for U = O.Ssf and 
Qi = O.Bep. /i is almost temperature independent below Tc, 
as is UeS- 



Schmitt-Rink (NSR)t^]. Ref. Q generalized the NSR 
approach to deal with a Feshbach resonance by including 
the contribution of 6-molecules outside the condensate 
(ignored in Eq. H1U|) ) in determining at Tc. 

In Ref. Ql, we extended the NSR approach to dis- 
cuss the superfluid phase below Tc in a uniform Fermi 
gas. Both A and /z are obtained by solving Eqs. © 
and (|10|l self-consistently, but now including the deple- 
tion of A through the presence of non-Bose-condensed b- 
molecules 22]. This procedure gives the simplest "renor- 
malization" of the MFA-BCS single-particle results in 
Eqs. ^ and which now involve values of fj, and A 
which include the effect of fluctuations around the MFA 
theory. The physics of this is most clearly shown in a 
functional integral treatment 23]. The gap equation © 
determines the MFA saddle point minimum, while the 
new number equation (|10|l describes the effect of Gaus- 
sian fluctuations around the MFA saddle point. The lat- 
ter gives rise to a renormalized value of the chemical po- 
tential ^ occurring in the MFA gap equation, which in 
turn leads to a renormalized value of A. It is important to 
note that this renormalized pairing approximation gen- 
erates response functions (^J which exhibit gapless Gold- 
stone modes, a crucial requirement of a conserving many- 
body approximation. Here we use the T — results to 
illustrate the behavior in a trapped Fermi gas 20] . 

In Fig. 1, we illustrate how Np and A'b vary as a 
function of the 6-moleculc threshold 2j/ for a trapped 
Fer mi g as at T = J(i]. The molecular condensate region 
Ii3,[l3ill3 occurs just below the two-body Feshbach res- 
onance at 2i/ — 0, where af^ is large and positive, leading 
to stable (long-lived) 6-molecules. In Fig. 1, the Fermi 
contribution includes the Cooper pairs. The latter con- 
tribution becomes less important for 2i/ < 2£p and is 
almost absent relative to the 6-molecules as soon as one 
passes through the Feshbach resonance (when v < 0, one 
has af > 0). A more detailed decomposition into Fermi 
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FIG. 2: Single-particle spectral density in a trapped gas of 
Fermi atoms at T = 0. The sharp peaks reflect the discrete 
BCS-Bogoliubov energies in a harmonic potential. 



excitations. Cooper-pairs and 6-molecules is given at Tc 
in Ref. 0, which may be more relevant to the region 
just below Tc studied in Refs. dill [13. 

In our pairing approximation, the effective attractive 
interaction Ucs in Eq. © plays a crucial role. One finds 
that Ucff smoothly increases in magnitude as v decreases, 
with nothing unusual happening at the two-body Fesh- 
bach resonance at 2z^ = Ofsee the inset in Fig. 1). 

As discussed in Ref. Q], the BCS-Bogoliubov quasi- 
particle spectral density is easily calculated from 



i Vlni[Gii(p,c^ + iO+)]. 

TT ^ ^ 



(14) 



The spectral density for a uniform gas at T = is shown 
in Fig. 8 in Ref. [3|. In the BEG region 2v < 0, we 
have fi < and as a result, /9f(^) has an excitation 



energy gap which starts at uj = y /i^ + |Ap for exci- 
tations with p = oH 0- In the BCS region 2v > Q 
(where ^ > 0), the quasiparticle gap starts at = |A|, 
from Fermi quasiparticles with momentum p = \/2m/i. 
In Fig. 2, we compare the density of states Pf(w) at 
V — e-p (BCS) with the result at the Feshbach resonance 
v = (BEC) in a trap. The sharp peaks correspond 
to the BCS-Bogoliubov excitation frequencies obtained 
from numerically solving(2Qj the Bogoliubov-de Gennes 
coupled equations, together with a self-consistent calcu- 
lation of and A(r) as described above. The sharp ex- 
citation edge in Fig. 2(a) reflects the BCS density of 
states (DOS) pf(w) ~ i?e[ ^ |A|^ "^' absent in 

the BEC region shown in Fig. 2(b). 

The threshold energy of the continuum spectrum in 
the density response function^ is equal to the minimum 
energy needed to break a bound pair of atoms and put 
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FIG. 3: Energy gap for atoms vs in a superfluid Fermi gas 
at T = trapped in a harmonic potential. 

them into the lowest energy states. This is the bind- 
ing energy Ebe of the dimer in the interacting Fermi 
gas. In the JILA experimentsfiol l2^ . a r.f. pulse of 
energy hv^.f is used to stimulate a transition to a lower 
Zeeman state. The two atoms will share an excess en- 
ergy AE = /iz^atom - hvr.f - Ebe, where /iJ^atom is the 
atom-atom transition frequency. If the dissociated pair 
of atoms are in the t (mp = ~9/2) and | (mp = —7/2) 
states involved in the molecular Bose condensate, the en- 
ergy of each atom (-Eatom) will have a threshold Eg given 

by £;atom > = \//^2 + |A|2 - in the BEC region 
{fi < 0), with the atoms having very small momentum. 
For A <C Ia*!, this gap Eg ~ |Ap/2|/x| is very small. 
The excess energy AE will be shared by the dissociated 
atoms, with -Eatom = Eg + In contrast, one has 

Eg = \A\ + n in the BCS region (/i > 0), with the atoms 
having a large momentum {p = y/2rnjl) comparable to 
the Fermi momentum (/i ^ ep ^ | A|). In Fig. 3, we plot 
Eg as a function of in a trapped Fermi gas at T — 0. 
Of course, these predictions assume that the dissociated 



atoms have thermalized with the rest of the atoms in the 
Fermi gas, so that the excitation spectrum is described 
by the equivalent of Eq. JHJ . 

By way of contrast, if the r.f. pulse dissociates the 
molecule into atoms in the spin state t (wf = —9/2) and 
another Zeeman state (mp = —5/2 in Ref. [i3)j only 
the t atom will be coupled to the underlying molecular 
condensate involving (|,|)-pairs. This atom will exhibit 
an energy gap Eg associated with the finite value of the 
order parameter |A|. The other atom will not have an 
energy gap. One sees that the energy of the pair of disso- 
ciated atoms can be different, assuming that a molecular 
condensate is only coupled to one of them. The energy 
and momentum spectrum of dissociated atoms should be 
a very direct probe of the pre sence (or absence) of a Bose 
condensate (see also Refs.[la,ll3 for other techniques to 
study single-particle excitations). 

In summary, we recall that in the original work on 
the BCS-BEC crossover 011, 0,0, the only bound states 
were Cooper pairs, whose existence is a many-body effect 
and does not require a two-body resonance. In this case, 
when tts > 0, the Cooper pairs are long-lived and form 
a molecular condensate. In contrast, in the presence of 
a Feshbach resonance, the condensate is formed of stable 
dimer molecules in the BEC limit, with a negligible con- 
tribution from Cooper pairs. We have studied the single- 
particle quasiparticle spectrum associated with a com- 
posite Bose order parameter arising from the formation of 
a BEC of Cooper pairs and Feshbach-induced molecules. 
We have emphasized that a simple static MFA pairing ap- 
proximation renormalized to include fluctuations 0,13 de- 
scribes both the BCS limit (where Cooper pairs dominate) 
and the BEC limit (where stable Feshbach molecules 
dominate). Both regions are described by a renormal- 
ized BCS-type single-particle Fermi spectrum. 
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